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ABSTRACT 

The  ability  to  rapidly  calculate  at-sensor  radiance  over  a  large  number  of  lines  of  sight  (LOSs)  is  critical  for 
hyperspectral  and  multispectral  scene  simulations  and  look-up  table  generation,  both  of  which  are 
increasingly  used  for  sensor  design,  performance  evaluation,  data  analysis,  and  software  and  systems 
evaluations.  We  have  demonstrated  a  new  radiation  transport  (RT)  capability  that  combines  an  efficient 
multiple-LOS  (MLOS)  multiple  scattering  (MS)  algorithm  with  a  broad-bandpass  correlated-/:  methodology 
called  ZURT-MS,  where  /dJRT  stands  for  corrclatcdWbascd  Ultra-fast  Radiative  Transfer.  The  MLOS 
capability  is  based  on  DISORT  and  exploits  the  existing  MODTRAN-D1SORT  interface.  ZURT-MS  is  a  new 
sensor-specific  fast  radiative  transfer  formalism  for  UV-visible  to  LWIR  wavelengths  that  is  derived  from 
MODTRAN’s  correlated-/:  parameters.  Scattering  parameters,  blackbody  and  solar  functions  are  cast  as  a 
few  sensor-specific  and  bandpass-specific  /'-dependent  source  terms  for  radiance  computations.  Preliminary 
transmittance  results  are  within  2%  of  MODTRAN  with  a  two-orders-of-magnitude  computational  savings. 
Preliminary  radiance  computations  in  the  visible  spectrum  are  within  a  few  percent  of  MODTRAN  results, 
but  with  orders  of  magnitude  speed  up  over  comparable  MODTRAN  runs.  This  new  RT  capability 
(embodied  in  two  software  packages:  ZURT-MS  and  MODTRAN-ZURT)  has  potential  applications  for 
remote  sensing  applications  such  as  hyperspectral  scene  simulation  and  look-up  table  generation  for 
atmospheric  compensation  analysis  as  well  as  target  acquisition  algorithms  for  near  earth  scenarios. 

Keywords:  /dJRT,  ZURT-MS,  Correlated-/',  DISORT,  Multiple  Scattering,  Look-up-tables,  MODTRAN, 
MODTRAN-ZURT,  Hyperspectral  and  Multispectral  Scene  Simulation 

1.  INTRODUCTION 

Look-up  tables  of  radiative  transfer  (RT)  quantities  such  as  bandpass  transmittance  and  radiance  are  used  in  a  variety  of 
remote  sensing  applications,  including  scene  simulation  and  atmospheric  compensation  of  spectral  imagery.  Since  the 
calculations  can  be  quite  time-consuming  if  performed  over  a  wide  range  of  atmospheric  conditions  and  viewing 
geometries,  highly  approximate  methods  are  often  used.  For  example,  the  Target  Acquisition  Weapons  Software 
(TAWS)  uses  the  2-stream  delta-Eddington  plane -parallel  atmosphere  multiple  scattering  (MS)  method. 12  This  method 
can  show  errors  of  30-50%  with  respect  to  the  “exact”  DISORT  (Discrete  Ordinates  Radiative  Transfer)  scattering 
algorithm3’4  and  displays  singularities  at  line-of-sight  (LOS)  and  solar  zenith  angles  of  90°;  i.e.,  for  horizontal  paths 
(HLOS)  and  the  sun  near  the  horizon.  For  such  broad  bandpass  applications,  and  other  applications  such  as  accurate 
hyperspectral  and  multispectral  scene  simulations,  there  is  both  the  need  and  opportunity  to  develop  new  fast  yet 
accurate  RT  methods  that  take  advantage  of  both  the  low  spectral  resolution  requirement  and  calculation  redundancies 
inherent  in  multiple  LOS  views  through  the  same  or  similar  model  atmospheres. 
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In  two  previous  SPIE  papers5'  6,  we  described  a  method,  dubbed  /<URT  (correlated-A-distribution  based  Ultra-fast 
Radiative  Transfer),  which  is  a  new  fast  and  fairly  accurate  thermal  radiance  RT  algorithm  designed  for  scene 
generation,  look-up  table  creation  and  target  acquisition  needs.  Here,  “fast”  is  defined  as  two  or  more  orders  of 
magnitude  more  rapid  than  1  cm"1  MODTRAN  calculations,  and  “accurate”  is  defined  as  an  absolute  transmission  error 
no  greater  than  2%  and  typical  radiance  errors  of  5-10%  when  compared  to  bandpass-averaged  1  cm_1-resolution 
MODTRAN  calculations.  In  addition,  /dJRT  automatically  accounts  for  the  sensor  bandpass  response  function,  thus 
eliminating  post-processing  of  higher-resolution  spectral  information.  The  second  SPIE  paper6  described  expanding 
kURT  ideas  into  computing  multiple-scattered  radiance  by  using  DISORT  within  the  /dJRT  framework,  dubbed  ZdJRT- 
MS.  This  paper  serves  as  an  update  on  both  the  previous  two.  Here  we  mainly  present  results,  leaving  out  details  of  the 
/dJRT- MS  formalism,  for  which  the  reader  is  referred  to  the  aforementioned  articles. 

The  spectral  range  of  /dJRT-MS  is  identical  to  that  of  MODTRAN;  i.e.,  from  the  infrared  (IR)  to  the  visible  (VIS). 
Preliminary  /dJRT-MS  computations  in  the  VIS-SWIR  indicate  that  RT  quantities,  particularly  radiance,  which  is 
dominated  by  MS  solar  radiance,  will  have  sufficient  accuracy  for  many  applications.  Of  course,  /dJRT-MS  will  not  be 
accurate  in  all  regions  of  the  spectrum,  especially  in  regions  of  low  transmittance.  Our  computations  in  the  visible  region 
indicate  that  the  kURT-MS  method  using  DISORT  predicts  total  radiance  for  HLOS  paths  with  low  solar  elevation 
angles  much  more  accurately  (2-6%  error)  than  the  2-stream  scattering  codes  (+30%  to  -55%  error). 

In  addition  to  /dJRT -MS,  which  can  run  independently  of  MODTRAN,  we  are  also  incorporating  /dJRT  concepts  into 
MODTRAN.  This  code,  which  we  call  MODTRAN-kURT,  will  be  fully  compatible  with  MODTRAN  and  can  run 
either  in  the  /dJRT  or  the  traditional  MODTRAN  mode.  M O D TRAN- k U R T  will  take  full  advantage  of  MODTRAN 
infrastructure  and  machinery,  including  I/O  and  refractive  geometry. 

Concurrent  with  /dJRT  development,  we  have  also  been  working  on  the  MODTRAN’s  DISORT  interface  to  take 
advantage  of  DISORT’s  inherent  multiple  line-of-sight  (MLOS)  capability,  namely,  the  capability  to  perform  RT 
calculations  for  many  lines  of  sight  using  a  single  call  to  DISORT  from  within  MODTRAN.  We  have  found  that  for  a 
computation  dominated  by  DISORT  (which  is  typical  for  a  calculation  covering  a  large  wavelength  region),  the  MLOS 
method  adds  only  15%  additional  time  for  an  additional  LOS  computation. 

2.  OVERVIEW  OF  THE  kURT-MS  METHOD 


2.1  Methodology  of  AURT-MS 

The  /dJRT  method  is  an  adaptation  of  the  traditional  correlated-/:  method  as  implemented  in  MODTRAN7  8' 9,  but  for 
“large”  bandpasses  and  with  the  instalment  sensor  relative  response  function  (RSR)  folded  into  the  RT  quantities.  The 
/dJRT  method  uses  pre-calculated  bandpass  specific  absorption  coefficient  values  (/(-values)  to  efficiently  compute 
bandpass  transmittances.  The  bandpass  specific  /(-values  are  created  off-line  and  are  the  results  of  “compacting”  a  large 
number  of  MODTRAN-generated  1  cm"1  /(-values  (called  fine-resolution  k- data  or,  simply,  Mb  for  /(-database,  typically 
17  per  cm"1)  to  a  mere  handful  of  /(-values  (typically  17  for  the  entire  band).  This  compact-/:  database  (cMb)  approach 
(see  the  left  side  of  Fig.  1)  reduces  the  number  of  monochromatic  calculations  for  a  broad  bandpass  from  thousands 
down  to  a  handful.  Note  that  MODTRAN  had  to  be  modified  to  output  the  Mb. 
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Fig.  1.  Major  Components  of /-URT-MS. 


More  on  the  Mb  and  cMb,  including  their  generation,  use  and  dependence  on  atmospheric  pressure  and  temperature  is 
described  in  the  earlier  SPIE  papers5  s.  There  we  also  describe  how  blackbody,  aerosol.  Rayleigh  and  cloud  scattering 
source  terms  (extinctions  coefficients,  phase  functions,  single  scattering  albedo)  are  correlated  with  cMb  and  with  the 
effective  transmittance  of  each  layer. 

As  described  in  the  aforementioned  references,  the  ATJRT  method  divides  the  atmosphere  into  important  absorbing 
species,  whose  concentrations  change  with  altitude,  and  the  rest  of  the  atmosphere  consisting  of  the  uniformly-mixed 
species  whose  relative  number  densities  do  not  change  with  altitude.  The  list  of  species  are  03,  C02,  H20  and  UMIX 
(uniformly  mixed),  the  last  of  which  includes  CO,  N2  continuum,  N20,  CH4,  NO,  02,  S02,  N02,  NH3  and  HN03. 
Treatment  of  H20  deals  with  H20  self-continuum  in  a  manner  different  from  H20  line-center,  line-tail  and  foreign- 
continuum  contributions.  Whereas  normally  optical  depth  is  proportional  to  absorber  amount,  for  H20  self-continuum 
the  optical  depth  scales  quadratically  with  concentration.  Thus,  we  have  five  molecular  “species”  in  ATJRT:  03.  C02, 
H20  without  self-continuum  contributions,  H20  self-continuum  and  UMIX. 


Assuming  a  lack  of  spectral  correlation  among  the  species,  the  transmittance  due  to  multiple  species  is  equal  to  the 
product  of  the  transmittances  of  the  individual  species.  With  17  cMb  terms  per  species,  this  leads  to  17s  terms  for  each 
layer  transmittance.  This  would  lead  to  175  c-k  terms  for  radiance  calculations  since  each  term  in  the  c -k  scheme  is 
treated  as  a  monochromatic  radiance  term.  Thus  there  is  a  need  to  reduce  the  number  of  c-k  terms  in  each  layer 
transmittance.  We  have  implemented  a  convolution-resampling  (C-R)  scheme  to  reduce  the  number  of  c-k  terms  for  this 
purpose  as  described  in  references5' 6.  The  method,  shown  in  Fig.  2,  begins  with  sorting  the  175  c-k  terms  and  concludes 
with  the  resampling  of  175  A:- values  to  just  17;  the  process  parallels  the  formation  of  cMb  from  the  Mb  data  shown  in 
Fig.  1.  In  practice,  instead  of  sorting  175  terms,  we  perform  the  C-R  operation  in  a  pair-wise  fashion,  which  involves 
performing  the  C-R  operation  on  4  xl72  terms,  a  factor  of  1228  smaller.  For  multiple  layers,  multiple  species 
transmittance  is  calculated  by  assuming  correlation  across  layers  of  the  optical  depth  Kj  L  C-r- 
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Fig.  2.  The  Convolution  and  Resampling  (C-R)  Scheme  to  Determine  Layer  Transmittance  of  the  Entire  Atmosphere.  The  reader 
is  referred  to  references5’ 6  for  further  detail,  including  explanation  of  mathematical  symbols. 

2.2  Anticipated  Use  of  AURT-MS  software  modules 

An  overview  of  the  anticipated  use  of  kURT  modules  is  given  in  Fig.  3.  Each  “frame”  corresponds  to  a  calculation 
process,  the  results  of  which  are  stored  for  later  use  by  the  inner  frame  processes.  An  outer  frame  denotes  a  process  that 
is  less  frequently  performed.  The  innermost  frame  consists  of  the  final  LOS  transmittance  or  radiance  calculation  that  is 
performed  repeatedly. 

The  outermost  process  is  the  one-time  assembly  of  the  full  high  spectral  resolution  /.-distribution  database,  covering  the 
entire  spectral  range  of  potential  interest.  The  next  process  represents  the  calculation  of  the  bandpass  ^-distribution  for 
each  atmospheric  layer  and  each  species.  This  process  may  be  broken  down  into  two  steps,  the  bandpass  convolution 
step  on  a  grid  of  pressures  and  temperatures,  and  for  the  case  of  a  specific  atmosphere,  interpolation  of  the  results  with 
respect  to  pressure  and  temperature.  For  a  standard  MODTRAN  model  atmosphere,  pre-interpolations  of  the  full  high- 
resolution  database  may  be  created,  avoiding  the  second  step.  The  third  process  involves  scaling  the  ^-distributions 
according  to  the  species  concentration.  In  addition,  for  the  calculation  of  radiance  we  implement  the  C-R  scheme.  In  the 
innermost  box,  the  correlated-A  approach  is  used  to  compute  path  radiances  and  transmittances. 

Once  fully  developed,  the  AURT-MS  package  will  be  a  suite  of  RT  routines  that  can  be  called  from  a  variety  of  software 
programs  such  as  TAWS  as  shown  in  Fig.  4.  Although  AURT-MS  is  being  written  in  FORTRAN  77  and  FORTRAN  90, 
wrappers  can  be  written  to  call  AURT-MS  routines  from  programs  in  any  language.  Currently,  we  are  developing  IDL 
(Interactive  Data  Language)  wrappers  for  AURT-MS. 
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Fig.  3.  Overview  of  the  kURT  calculation  scheme. 


Fig.  4.  Schematic  of  the  HJRT-MS  interface  architecture.  The  calling  program  will  access  the  /rURT-MS  library  through  a 
collection  of  wrapper  routines. 


3.  RADIANCE  AND  TRANSMITTANCE  RESULTS 

We  have  used  the  kURT  method  to  compute  transmittances  and  radiances  in  the  3-5  pm,  8-12  pm  and  visible  wavelength 
regions  for  down-looking  and  horizontal  paths.  The  3-5  pm  region  has  been  the  most  difficult  region  to  model  because  of 
the  strong  CCR  absorption  band  from  4. 2-4. 5  pm  and  strong  FTO  absorption  from  4. 5-5. 2  pm.  kURT-MS,  at  this  point  of 
development,  gives  poorest  results  when  thermal  radiance  is  both  small  and  the  dominant  contributor  to  radiance  ( e.g .,  in 
the  3-5  pm  region)  and  when  the  transmittance  is  small  but  not  negligible.  The  best  results  are  when  solar  scattering 
dominates,  as  in  the  VIS-SWIR  region.  There  is  always  excellent  agreement  with  ATJRT  and  MODTRAN  transmittance 
values. 


3.1  Down-looking,  10-0  km  path,  radiance  and  transmittance  results  in  the  3-5  pm  region 

Tables  1  and  2  present  RT  results  for  10  km  down-looking  paths  for  the  six  MODTRAN  model  atmospheres.  The  model 
atmospheres  range  from  the  wettest  tropical  to  the  driest  Sub-Arctic  Winter.  The  vertical  H20  amounts  are:  Tropical 
(model  1)  5120  atm-cm,  Mid-Latitude  Summer  (model  2)  3636  atm-cm,  Sub-Arctic  Summer  (model  4)  2590  atm- 
cm,  US  Standard  (model  6)  1762  atm-cm,  Mid-Latitude  Winter  (model  3)  1060  atm-cm,  and  Sub-Arctic  Winter 
(model  5)  518  atm-cm.  The  results  are  fairly  good  for  the  C02  absorption  band  (Table  2).  For  unknown  reasons  the 
radiance  results  tend  to  get  worse  in  terms  of  percentage  error  as  the  F^O  amount  in  the  atmosphere  decreases  (Table  1). 


Table  1.  5.2-4. 5  pm  Bandpass  (1925-2220  cm  ')  Thermal  Radiances  and  Transmittances  for  the  10-0  km  Vertical  Paths  for 
the  Six  MODTRAN  Model  Atmospheres.  For  all  LOSs,  C02  concentration  is  365  ppmv. 
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TROP 

MLS 

SAS 
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MLW 

SAW 

5.2-4.5  pm 

Model  1 

Model  2 

Model  4 

Model  6 

Model  3 

Model  5 

MOD  Rad 

6.227E-05 

5.1 10E-05 

3.376E-05 

2.909E-05 

1.658E-05 

8.122E-06 

(w/sr/cm2),  T 

0.257 

0.317 

0.387 

0.444 

0.522 

0.613 

AURTT 

0.274 

0.349 

0.426 

0.489 

0.566 

0.645 

T  Error 

0.017 

0.032 

0.039 

0.045 

0.044 

0.032 

AURTR 

5.79E-05 

4.56E-05 

2.93E-05 

2.45E-05 

1.34E-05 

6.33E-06 

R  Error  (%) 

-6.97 

-10.84 

-13.36 

-15.94 

-19.43 

-22.02 

Table  2.  4. 5-4. 2  pm  (2220-2380  cm  *)  Bandpass  Thennal  Radiances  and  Transmittances. 
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TROP 

MLS 

SAS 

USS 

MLW 

SAW 

4. 5-4. 2  pm 

Model  1 

Model  2 

Model  4 

Model  6 

Model  3 

Model  5 

MOD  Rad 

4.027E-06 

3.586E-06 

2.320E-06 

2.005E-06 

1.485E-06 

9.874E-07 

(w/sr/cm2),  T 

0.009 

0.010 

0.012 

0.011 

0.012 

0.013 

AURTT 

0.024 

0.026 

0.028 

0.029 

0.030 

0.032 

T  Error 

0.015 

0.016 

0.016 

0.018 

0.018 

0.019 

AURTR 

4.24E-06 

3.74E-06 

2.36E-06 

2.08E-06 

1.51E-06 

9.61E-07 

R  Error  (%) 

5.28 

4.22 

1.77 

3.86 

1.68 

-2.68 

3.2  Down-looking,  10-0  km  path,  radiance  and  transmittance  results  in  the  8-12  pm  region 

Tables  3  and  4  present  RT  results  for  10  km  down-looking  paths  for  the  six  MODTRAN  model  atmospheres.  FI  ere  the 
both  transmittance  and  radiance  results  are  satisfactory,  although,  again,  the  percentage  radiance  errors  increase  as  the 
FLO  content  of  the  atmosphere  decreases. 


Table  3.  10.664-8  pm  (940-1250  cm  ')  Bandpass  Thermal  Radiances  for  the  10-0  km  Vertical  Paths  for  the  Six  MODTRAN 
Model  Atmospheres. 


10-0  km 

TROP 

MLS 

SAS 

USS 

MLW 

SAW 

10.664-8  pm 

Model  1 

Model  2 

Model  4 

Model  6 

Model  3 

Model  5 

MOD  Rad 

9.138E-04 

6.258E-04 

3.792E-04 

2.677E-04 

1.549E-04 

7.931E-05 

(w/sr/cm2),  T 

0.544 

0.649 

0.736 

0.803 

0.847 

0.892 

AURTT 

0.551 

0.659 

0.745 

0.808 

0.854 

0.897 

T  Error 

0.007 

0.010 

0.009 

0.005 

0.007 

0.005 

AURT  R 

9.02E-04 

6.11E-04 

3.68E-04 

2.62E-04 

1.49E-04 

7.61E-05 

R  Error  (%) 

-1.25 

-2.42 

-2.88 

-2.02 

-4.00 

-4.06 

Table  4.  8-7.377  pm  (1250-1355  cm  ')  Bandpass  Thermal  Radiances  and  Transmittances  for  the  10-0  Vertical  Paths  for  the 
Six  MODTRAN  Model  Atmospheres. 


10-0  km 

TROP 

MLS 

SAS 

uss 

MLW 

SAW 

8-7.377  pm 

Model  1 

Model  2 

Model  4 

Model  6 

Model  3 

Model  5 

MOD  Rad 

2.628E-04 

2.436E-04 

1.953E-04 

1.797E-04 

1.388E-04 

9.360E-05 

(w/sr/cm2),  T 

0.029 

0.049 

0.079 

0.110 

0.167 

0.257 

ATJRTT 

0.029 

0.055 

0.093 

0.135 

0.204 

0.298 

T  Error 

0.000 

0.006 

0.014 

0.025 

0.037 

0.041 

ATJRTR 

2.75E-04 

2.52E-04 

1.99E-04 

1.81E-04 

1.35E-04 

8.92E-05 

R  Error  (%) 

4.62 

3.28 

1.65 

0.60 

-2.73 

-4.75 

3.3  10  km  near-HLOS  path  RT  results  in  the  0.4-0.9  pm  (VIS-SWIR)  bandpass 

Table  5  shows  transmittance  RT  results  for  a  10  km  near-HLOS  path  (from  0.1  km  to  0.1  km)  using  the  US  Standard 
MODTRAN  atmosphere  and  two  aerosols  models.  The  AXJRT  results  for  the  wavelength  region  0.4-0. 9  pm  are  excellent. 
The  two  regions  correspond  to  11,000-15,000  cm'1  and  15000-25,005  cm"1. 

Table  5.  Single-  and  Two-Region  Transmittance  Values  for  a  10  km  near  HLOS  path  using  the  US  Standard  MODTRAN 
model  atmosphere  and  two  aerosol  models. 


Aerosol,  PATH  LENGTH  and  O.D. 

2-region, 

AURT 

1-region, 

AURT 

MODTRAN 

Rural,  10  km,  OD=2 

0.1388  (+0.003) 

0.1447  (+0.008) 

0.1359 

Rural,  10  km,  OD=3 

0.0597 

0.0625  (+0.003) 

0.0588 

Urban,  10  km,  OD=2 

0.1364 

0.1420  (+0.008) 

0.1336 

Urban,  10  km,  OD=3 

0.0576 

0.0603  (+0.003) 

0.0567 

As  with  the  transmittance  calculations,  we  performed  AURT  radiance  calculations  using  single-region  and  two-region 
methods.  The  two-region  results  are  significantly  better  than  the  one-region  results.  A  typical  set  of  2-region  results  is 
presented  in  Fig.  5.  Here  the  solar  zenith  angle  is  80°.  There  are  12  pairs  of  calculations.  Each  pair  consists  of  a  2-stream 
MODTRAN1011  calculation  (purple  and  negative  bars)  and  AURT  DISORT  calculation  (dark  red  and  errors  are  all 
positive).  The  first  six  pairs  correspond  to  rural  aerosol,  optical  depth  of  2  and  3,  and  the  next  six  pairs  correspond  to 
urban  aerosols  of  optical  depth  2  and  3.  Each  grouping  of  three  pairs  corresponds  to  azimuths  of  180°  (staring  away  from 
the  sun),  90°  and  0°  (staring  at  the  sun).  As  expected,  when  looking  towards  the  sun  2-stream  errors  are  at  their  peak.  In 
general  /TJRT-DISORT  errors  versus  MODTRAN  with  8-stream  DISORT  are  less  than  6%,  whereas  2-stream  errors  are 
as  high  as  55%. 


Fig.  5.  Isaacs  2-stream  and  kURT-MS  0. 4-0.9  pm  radiance  comparisons  to  MODTRAN  DISORT  (“ground  truth”)  for  the  ARL 
scenario.  Results  are  shown  for  rural  (cases  1-6)  and  urban  (7-12)  aerosols  of  optical  depth  2  (cases  1-3,  7-9)  and  3  (4-6, 
10-12).  Each  triple  compares  results  for  backward,  broadside  and  forward  scattering.  The  solar  zenith  angle  is  80°. 


3.4  Down-looking  visible-MWIR  radiance 

We  carried  out  several  off-nadir  VIS-MWIR  RT  calculations  (Table  6)  using  the  mid-lat  summer  MODTRAN 
atmosphere  with  H20  scaled  to  a  fairly  high  value  (2.5  cm,  Scenario  1)  and  the  winter  atmosphere  with  H20  scaled  to  a 
fairly  low  value  (0.4  cm,  Scenario  2).  The  aerosol  was  rural  with  visibilities  of  23  and  5  km.  Ground  reflectance  was 
10%  and  the  C02  was  365  ppmv.  The  solar  and  observer  elevation  angles  are  both  45°  above  the  horizon,  and  the  sun  is 
located  in  front  of  the  observer  (MODTRAN  azimuth  is  0°).  We  considered  several  generic  (not  top-hat)  response 
functions.  Results  from  8-stream  MODTRAN-/. U R T  calculations  are  shown  below.  Only  the  highest  errors  are  shown 
within  parenthesis,  and  are  omitted  when  they  are  negligible.  The  highest  radiance  error  is  +4%,  but  typically,  much 
lower.  For  these  wavelengths  and  scenarios,  it  is  the  solar  scattering  that  constitutes  almost  all  of  the  total  radiance. 

We  also  performed  timing  analyses  on  these  types  of  calculations.  We  found  that  for  a  4000  cm"1  AURT  calculation, 
AURT  was  faster  by  a  factor  15  over  15  cm"1  MODTRAN  band  model  calculation  and  about  200  times  faster  than  1  cm"1 
MODTRAN  band  model  computations.  These  timing  ratios  are  consistent  with  the  number  of  spectral  calculations 
performed  by  the  respective  methods. 

Table  6.  VIS-MWIR  Nadir  Calculations. 


Scenario  1  Scenario  2 


FILTER  (pm) 
Effective  Width 
Visibility 

Transmittance 

Total  Rad 
(pW  cm-2  sr1) 

Transmittance 

Total  Rad 
(pW  cm-2  sr+) 

A  0.438  -  0.536 
2092  cm'1 

23  km 

AURT 

0.4643  (+0.0013) 

406 

0.4865 

393 

M5 

0.4630 

403 

0.4853 

390 

5  km 

AURT 

0.1191 

572 

0.1249 

564  (+0.9%) 

M5 

0.1187 

567 

0.1246 

559 

F  0.895  -  0.985 
611.6  cm-1 

23  km 

AURT 

0.3653 

38.8  (+0.5%) 

0.6265  (+0.0069) 

69.6  (+2.4%) 

M5 

0.3647 

39.0 

0.6196 

68.0 

5  km 

AURT 

0.2021 

56.0  (+2.0%) 

0.3435  (+0.0053) 

104  (+4.0%) 

M5 

0.1994 

54.9 

0.3382 

100 

D  0.730-0.889 
1101  cm-1 

23  km 

A 

0.6747  (+0.0022) 

184  (+0.5%) 

0.7209 

(+0.0005) 

189 

M  5 

0.6725 

183 

0.7204 

189 

5  km 

A 

0.3199  (+0.0023) 

282  (+0.7%) 

0.3401 

(+0.0006) 

298 

M  5 

0.3176 

280 

0.3395 

297 

H  1.352-1.416 
139.5  cm  1 

23  km 

A  UR  1 

0.0000 

0.274  (-2.5%) 

0.0427(+0.0012) 

0.540  (-3.0%) 

M  5 

0.0000 

0.281 

0.0415 

0.557 

5  km 

/.l  R I 

0.0000 

0.274  (-2.5%) 

0.0303  (-0.0013) 

0.681  (-2.1%) 

M  5 

0.0000 

0.281 

0.0316 

0.696 

J  3.30-4.26 
381.9  cm  1 

23  km 

AIR  1 

0.7490  (+0.0097) 

23.8  (+0.8%) 

0.8279(+0.0003) 

13.9  (+3.0%) 

M  5 

0.7393 

23.6 

0.8282 

13.5 

5  km 

AURT 

0.6603  (+0.0082) 

24.2  (+1.7%) 

0.7272  (+0.0002) 

14.5 

M  5 

0.6521 

23.8 

0.7270 

14.5 

4.  MODTRAN  SPEED-UP  USING  DISORT  MLOS  STRATEGY 


The  DISORT  code  outputs  information  that  can  be  used  to  compute  radiances  for  multiple  lines-of-sights  (MLOSs)  with 
a  single  subroutine  call.  The  one  caveat  is  that  either  the  sun  and  sensor  latitudes  and  longitudes  or  the  sun  and  target 
latitudes  and  longitudes  must  be  fixed  for  the  entire  set  of  calculations.  However,  the  MODTRAN-DISORT  interface 
has  not  until  now  been  programmed  to  take  advantage  of  this  capability.  The  basic  concept  within  DISORT  is  reasonably 
straightforward.  Given  a  stratified  plane-parallel  atmosphere  and  a  fixed  solar  geometry,  the  discrete  ordinate  MS 
algorithm  DISORT  solves  the  monochromatic  radiation  transport  equation  (RTE).  Output  includes  mth  azimuth  moment 
radiant  intensities,  I  m(  r,  ±ju),  at  user-specified  nadir  optical  depths,  r.  within  the  atmosphere  and  at  user-specified  up- 
look  ( +ju )  and/or  down-look  (-//)  view  angles  (/j  >  0  is  the  cosine  of  the  zenith  or  nadir  view  angle,  respectively).  The 
final  azimuth  (^-dependent  radiant  intensities  are  computed  from  a  truncated  Fourier  cosine  series: 

2M-1 

Im(T,ju,</>)=  ^  (r ,  //)  cos  m(<f>0  -  <f)  (1) 

m= 0 

where  (f>o  is  the  solar  azimuth  angle  and  2 M  is  the  number  of  azimuth  moments. 

The  MODTRAN-DISORT  MLOS  option,  when  fully  implemented,  will  greatly  accelerate  the  generation  of  MLOS 
look-up  tables  by  exploiting  the  full  set  of  results  acquired  from  a  single  run  of  DISORT.  For  single-LOS  MODTRAN 
runs  that  include  DISORT  MS  and  a  large  number  of  spectral  grid  points,  the  DISORT  module  execution  time 
completely  dominates  (>99%)  the  overall  run  time.  However,  the  DISORT  MS  results  for  an  additional  view  angle  can 
be  obtained  with  about  a  15%  increase  in  the  baseline  execution  time.  Thus,  the  MLOS  option  enables  the  users  to  run 
~8  view  angles  in  the  same  time  that  the  single  LOS  MODTRAN  would  perform  2  view  angle  calculations. 
Furthermore,  varying  path  azimuth  angle,  sensor  altitude  and/or  target  altitude  adds  less  than  1%  to  the  baseline 
execution  time.  Thus,  a  grid  of  8  view  angles,  8  path  azimuth  angles,  8  sensor  altitudes  and  8  target  altitudes,  a  total  of 
more  than  4000  LOSs,  will  run  in  about  the  same  time  as  3  single-LOS  MODTRAN  runs. 

This  capability  is  also  ideally  suited  for  3D  radiance  scene  simulations,  which  typically  involve  radiance  and/or  contrast 
signature  imagery  over  a  short  time  duration.  Thus,  sun  location  relative  to  the  Earth  is  essentially  constant,  and 
radiation  transport  solutions  may  be  required  for  thousands  of  LOSs.  In  the  past,  accurate  incorporation  of  the  MS  signal 
has  been  thought  to  be  computationally  prohibitive.  However,  if  the  myriad  of  views  are  extracted  from  a  single  MS 
calculation,  the  computation  of  high-fidelity  background  scattered  radiances  for  each  imaged  pixel  is  an  achievable  goal. 

5.  CONCLUSION 

An  instrument  response  function-specific  correlated-^  approach,  dubbed  /dJRT-MS.  is  being  developed  to  rapidly 
compute  bandpass  LOS  transmittances  and  radiances  from  the  visible  to  the  IR.  In  its  frill  implementation,  /dJRT -MS 
will  be  a  complete  RT  code  that  includes  scattering  due  to  aerosols  and  clouds.  /dJRT -MS  software  will  be  a  complete 
library  of  RT  routines  that  includes  all  significant  radiance  components;  and  the  individual  routines  of  the  library  can  be 
called  on  an  a-la-carte  basis  for  a  variety  of  applications  that  require  forward  RT  modeling.  The  goal  is  a  hundred-  to 
thousand-fold  speed-up  with  2%  accuracy  in  transmittance  and  10%  accuracy  in  radiance  compared  to  a  1  cm-1 
MODRTRAN  run.  Preliminary  results  on  several  wide  bandpasses  and  narrower  visible-MWIR  instrument  response 
functions  indicate  that  both  accuracy  and  speed  goals  are  within  reach.  We  are  also  incorporating  the  HJRT 
methodology  into  a  version  of  MODTRAN,  called  MODTRAN-AURT ,  which  will  run  from  a  MODTRAN  interface 
either  in  the  AURT  mode  or  in  the  usual  MODTRAN  mode. 

Additionally,  MODTRAN’s  DISORT  interface  is  being  upgraded  to  take  advantage  of  DISORT’s  inherent  multiple  line- 
of-sight  (MLOS)  capability,  namely,  the  capability  to  perform  RT  calculations  for  many  lines  of  sight  using  a  single  call 
to  DISORT  from  within  MODTRAN.  This  will  result  in  very  substantial  savings  in  computation  time  for  scene 
modeling  and  radiance  look-up  table  generation. 
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